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^ ' Abstract 

We construct boson stars in global Anti de Sitter (AdS) space and study their stability. 

Linear perturbation results suggest that the ground state along with the first three 

excited state boson stars are stable. We evolve some of these solutions and study their 
\Q • nonlinear stability in light of recent work [T] arguing that a weakly turbulent instability 

drives scalar perturbations of AdS to black hole formation. However evolutions suggest 
. that boson stars are nonlinearly stable and immune to the instability for sufficiently 

Q ! small perturbation. Furthermore, these studies find other families of initial data which 

cn I similarly avoid the instability for sufficiently weak parameters. Heuristically, we argue 

that initial data families with widely distributed mass-energy distort the spacetime 

sufficiently to oppose the coherent amplification favored by the instability. From the 
/S . dual CFT perspective our findings suggest that there exist families of rather generic 

c^ ! initial conditions in strongly coupled CFT (with large number of degrees of freedom) 

that do not thermalize in the infinite future. 
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1 Introduction 

Understanding the gravitational behavior of spacetimes which asymptotically 
behave as anti deSitter (AdS) has, since holography [2], attracted significant 
interest. The AdS/CFT correspondence conjecture implies that such under- 
standing is critically important for a plethora of phenomena described by field 
theories. Remarkably however, relatively little is known about dynamical 
scenarios, especially in comparison to spacetimes which are asymptotically 
fiat (AF) or asymptotically deSitter (dS). An important reason behind this 
difference is that the boundary of AdS is in causal contact with the interior 
of the spacetime, in stark contrast to the boundaries of AF and asymp- 
totically dS spacetimes. That the boundary affects the interior prevents a 
straightforward extension of standard singularity theorems to asymptotically 
AdS (aAdS) spacetimes [3]. 

Significant advances have recently been achieved via different approaches, 
including strictly analytic IH El [6l [7], perturbative [H |9] and numerical [H 
[T0| [TTl [T2| [T3| [HI [15] (to name a few representative) efforts. Particularly 
intriguing is the evidence first presented in [1] that pure AdS is unstable to 
scalar collapse to black hole (BH), regardless of how small a scalar perturba- 
tion is considered. The effect of the AdS boundary allows for the reflection 
of scalar pulses. Hence, a weak scalar pulse bounces off the boundary at 
inflnity, returning to concentrate again at the origin. The weakly turbulent 
instability of [1] results in the sharpening of the pulse so that, eventually, it 
achieves sufficient concentration to form a black hole. The black hole critical 
behavior discovered by Choptuik [16] in AF spacetimes appears here repeat- 
edly. In particular, after each reflection off the AdS boundary, there is yet 
another threshold for prompt BH formation (before the next bounce). 

To explain this behavior, a study of the normal modes for scalar [1] and 
tensor [01 perturbations revealed a non-linear mode coupling at third order, 
shifting energy to higher frequencies. Consequently, this shift in frequency 
would eventually lead to the formation of a black hole and so generic instabil- 
ity of AdS was conjectured. Interestingly though, further studies suggested 
the existence of non- linearly stable solutions [HI [131 [151 [II]- These obser- 
vations hint that a straightforward application of the resonance picture can 
only partially capture the dynamical behavior. 

We consider the dynamics of boson stars in AdS, motivated towards bet- 
ter understanding of this instability and its implications for holographic sce- 
narios. Horizon formation in gravitational collapse from the holographic 
perspective implies thermalization of the dual conformal theory. When the 
gravitational problem is formulated in global d + 1 space-time dimensional 
asymptotically AdSd+i, the dynamics of the dual d space-time dimensional 
CFTd occurs on a S'^~^ sphere. Furthermore, the holographic duality arises 
as a correspondence between a String Theory and a CFT — it can be trun- 



cated to a gravity/CFT duality when a CFT has a large number of strongly 
interacting degrees of freedom!^. A natural expectation from the field theory 
perspective is that a large number of strongly interacting degrees of freedom 
in a finite volume thermalize from generic initial condition; thus, one does 
expect no-threshold BH formation as advocated in [1]. Early work on bo- 
son stars in AdS [18] demonstrated that such stationary configurations are 
linearly-stable. Thus, initializing a CFT state as a dual to an AdS boson star 
might result in a slow thermalization of the latter. In this paper we show 
that, rather, the boson stars in AdS are non-linearly stable. This non-linear 
stability is unaffected by boson star perturbations, as long as these pertur- 
bations are sufficiently small. Even more surprising, we find that families of 
initial conditions with widely distributed mass-energy are non-linearly stable 
for sufficiently small mass as well. This suggests that non-linear stability is 
not a feature of states carrying a global charge (as dual to boson stars). As a 
result, it appears that there exist large sets of initial configurations in CFT 
which never thermalize in their evolution. 

The paper is organized as follows. In section [2] we setup the gravitational 
dual of the generic, spatially-isotropic CFT^, initial condition specified by a 
pair of dimension-3 operators with a global U{1) symmetry: an Einstein grav- 
ity with a negative cosmological constant and a massless, minimally coupled 
bulk complex scalar field. In section |3] we construct stationary configurations 
of the coupled scalar-gravity system, charged under the global f/(l) symme- 
try — boson stars. We show that boson stars are stable under linearized 
fiuctuations. In section H] we present results of fully-nonlinear simulations of 
genuine boson stars, perturbed boson stars, and f/(l)-neutral initial configu- 
rations with widely distributed, bulk mass-energy. We conclude in section \5\ 
and outline future directions. 



2 Effective action and equations of motion 



In this section, following [13], we review the formulation of the problem of 
the gravitational collapse of a complex scalar in asymptotically anti-de- Sitter 
space-time. We focus on asymptotically AdS4^ collapse, dual to CFT3 (we 
choose d = 3 in the notation of [13]). 

The effective four- dimensional action is given bjo 

"^4 = TTT?^ / ^'^v^ (i?4 + 6 - 2d,(Pdy*) , (2.1) 

167rG4 Jj^^ 



^This can be quantified as a large central charge for even d, or more generically, the 
large number of excited degrees of freedom at thermal equilibrium. 
^We set the radius of AdS to one. 



where = 0i + i 02 is a complex scalar field and 

Mi = dM-ixX, dM3 = RtxS\ l={xe[0,^]}. (2.2) 
Adopting the line element as in [19], 



^ (-Ae-^'dt' + ^' 
cos"^ X \ A 



ds^ = -— ^ ( -Ae-^^dt^ + -r + sin^ x dnl] , (2.3) 



where rffig is the metric of unit radius S*^, and A{x,t) and S{x,t) are scalar 
functions describing the metric. Rescaling the matter fields as in 



(2.4) 
n^ = T^^' (2-5) 

A COS^ X 

$. = ^, (2.6) 

cosx 

we find the following equations of motion (we drop the caret from here for- 
ward) 



COS^ X ' 


A cos^ X 
dx(t)i 


cosx ' 



--Ae-^Iii 



^'•■, 



cos X 

■ If sin^ X , c ^ 

n. =—^ Ae~^<l>i 

sm X \ cos X 



(2.7) 



l + 2sin2x ,, ,, . F. . f ^f 



5 „ A i I Tt2 



A ^ = (1 - A) - sin X COS*' xA — ^ + 11: 

sm X cos X \ COS"' x 



5 ^ = — sin X cos^ X I — ^ h n,^ 



cos^ X 



(2.8) 



together with one constraint equation 

At + 2 sin X cos^ A^e'^ ($^0,) = , (2.9) 

where a sum over i = {1,2} is implied. We are interested in studying the 
solution to (I2.7p -( l2l9|) subject to the boundary conditions: 
■ Regularity at the origin implies these quantities behave as 

A{t,x)=l + 0{x^), (2.10) 

6{t,x)=6o{t) + 0{x^)- 



at the outer boundary x = 7t/2 we introduce p = 7r/2 — x so that we have 



A{t,p)=l-M^^^ + Oip'] 
cosp 

6{t,p)=0 + O{p'). 



(2.11) 



The asymptotic behavior fl2.1ip determines the boundary CFT observ- 
ables: the expectation values of the stress-energy tensor T^i, and the opera- 



tors O3 , dual to (f)i, 



87rG,{Tu) = M, 
167rG.+i(O?) = 12 0«(t) 



(Ta,) = ^-f (Tu) , 



(2.12) 



where gap is a metric on a round S^. Additionally note that the conserved 
t/(l) charge is given by 



Q = 87r 



n/2 



dx sin^x cos^x (03(0, x)0i(O, x) - ni(O,x)02(O, x)) , (2.13) 



and that since dtQ = 0, the integral in ( 12.13^ can be evaluated at t = 0. 

The constraint (12. 9 p implies that M in (12. lip is time-independent, ensur- 
ing energy conservation 

dt{Tu) = 0. (2.14) 

It is convenient to introduce the mass aspect function Ai{t,x) as 



COS'^ X 



A{t,x) = 1-M{t,x 

smx 

Following (12. 8p we find 

Ai{t,x) = / dz tan^ 2; cos zA{t,z 
Jo 

Comparing (I2.16p and (12. lip we see that 



(2,15) 



^ + n^(,.) 



COS^ z 



(2.16) 



M = M{t,x) 



(2.17) 



3 Boson stars in AdS4 

There is an interesting class of stationary, perturbatively stable, fully-nonlinear 
solutions to f l2.7p -( l2rTT]) with nonzero Q, referred to as boson stars [T8l [20]. 
Such solutions are characterized by a discrete integer n = 0, 1, ■ ■ ■ , denot- 
ing the number of nodes of the complex scalar radial wave-function, and 
a continuous value of the global charge Q. In this section we discuss the 
numerical construction of such solutions, their perturbative properties, lin- 
earized stability, and their relation (for small Q) to linearized AdS^ massless, 
minimally coupled, scalar modes — the oscillons. Oscillons were reviewed 
in detail in [13]. These stationary solutions are uniquely characterized by an 
excitation level j = {0, 1, ■ ■ ■ }: 



Cj {x) = dj cos^ X 2F1 ( -j, 3 + j; -; sin^ x j , 



(3.1) 



where w^^^ is an oscillon frequency and dj is a constant enforcing their or- 
thonormality, 

/■7r/2 

/ dx ei{x)ej{x) tan^x = 6ij . (3.2) 

3.1 Stationary boson stars 

Assuming a stationary solution in which the complex field varies harmonically 



'f^i^) Ju^t 



(j)i(x,t)+i(j)2(x,t) = \ e"^\ A(t,x) = a(x) , 6(t,x)=d(x) 

cos"' a; 



(3.3) 
we find a system of ODEs from (12. 7p and (12.81) 

=4>" + ( + --d') 4>' + u'e^'a-' , 

\cosxsmx a J 

=d' + sin X cos X a-2 {{^'fa^ + cp^u^e^'') , (3.4) 

=a' + ^'^^ (1 -a) + sinxcosx a~^ ((^O^^^ + ^^^^e^"*) . 

cos X sin X 

The charge and the mass determined by these solutions are given by 



^ „ ;^/2 , wsin2x0(x)2e'^(^) 
Q =8n I dx ^—^ 







aix) cos'x 



^/2 sin^ X 



(3.5) 



^= I ^^ ..ZZZ2^ (a(x)2(0'(x))2 + e"'^^^uj^(t>{xy) 
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Q al^xjcos^x 



Physically relevant solution to equations 13.41 must satisfy: 

■ At the origin of AdS, i.e., x — )■ 0: 

a=i-l(pQ^)Ve2'^o x2 + 0(a:^), 
o 

d =d^^ - ^ip^,)We''o x' + Oix') , (3.6) 

o 

Note that besides u, the general solution is characterized by 

K,4}. (3.7) 

■ and asymptotically (at the AdS boundary, i.e., p -^ 0): 

a=l + alp' + 0{p') , 

d =\{l^,fp' + {^^{V\f - \{vlfuj') P' + 0{p') , ^3 g^ 

=pIp' + (^IpI - ^pI-^') p' + C?(/) • 

Note that besides u, the general solution is characterized by 

{p14}. (3.9) 

From fl3.7p and fl3.9p we have precisely the correct number of coefficients 
to find an isolated solutioiio for a given co. 

3.2 Small charge boson stars as oscillons 

In this section we discuss analytic results for the spectrum of boson stars 
perturbatively in the amplitude (p. To this end we introduce 

= A 01 + 0{X') , a = 1 + A^ a2 + 0{X') , 

e-^ = 1 + A' ^2 + C(A^) , cu = cuo + U2\^ + O(A^) , 

where A is an expansion parameter. Substituting fl3.10p into (13. 4p we find 
that equation for 0i decouples to leading order 

=01 + —^ 0'i + u:l 01 . (3.11) 

sm X cos X 



■^As we discuss in section 13.21 solutions of a fixed charge are labeled by an integer, 
specifying the "level" of a boson star. 
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Normalizing 0i as 

= 1, (3.12) 

the general solution of (13.111) . subject to the boundary conditions (13. 6p and 
(13.81) . is given by 

0S'^ = j- ej{x) , 4') = w^'^ = 3 + 2j , (3.13) 

Uj 

where the integer j = 0, 1,2--- parameterize the 'excitation level' of the 
boson star. Notice that the asymptotic expansion for (in eqn 13. 6p implies 
Pq = X . Given (I3.13p . we find from (13. 5p 

""-^ 8(j + l)(j + 2)+^^^)' ^-' 2(, + l)(, + 2)+^^^)' 



(3.14) 



Note from (13.140 that for a fixed charge Q, excited levels of boson stars are 
more massive. 

It is straightforward to compute the leading-order background warp fac- 
tors {02,(^2}, as well as subleading frequency correction U2. In what follows 
we present explicit expressions for the first four levels of a boson star: 
■ j = level, 

df^ = lcos'^x, (3.16) 



(0) 63 



2 

~32 
j = 1 level, 



(3.17) 



^(1) ^ 25 cos X /I in(4a;)(2cos(4x) + 1) - 3x ) , (3.18) 



72 sin X V 4 



7(1) 5 



4^^ = _ cos^ X (32 cos^ X - 40 cos^ x + 15) , (3.19) 
18 

a;« = -^; (3.20) 

2 864 ' ^ ^ 



j 


"7/' + 0(A2) 


i + o{\') 


i + o{x^) 


|| + 0(A2) 





63 
32 


3 

2 


1 


Q 


1 


3025 
864 


35 

18 


5 
3 


-i- 


2 


89327 
17280 


203 
90 


7 
3 


-i- 


3 


154143 
22400 


873 
350 


-3 


160 " 



Table 1: Condensate values (13. 6 p and ( 13. 8 p of boson stars, perturbatively in 
in Pq = A. 



j = 2 level, 



.(2) 



49 cos^ X / 1 



144 sin X V8 



sin(8x)(2 cos(4x) + 1) — 3x 



(3.21) 
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dr = — cos^ X (960 cos^ X - 2240 cos^ x + 1904 cos^ x - 700 cos^ x + 105) 
"^ 90 ^ ' 



w 



(2) 



89327 _ 
"17280' 



(3.22) 
(3.23) 



j = 3 level. 



l'S^ 81 COS X I I , ^ , , X , X 

a^'^^ = — — I - sin(8x)(2 cos(8x) + 2 cos(4x) + 1) - 5x 

TlUU olll Jb \ O 



(3.24) 



4^^ = cos^ x(28672 cos^^ x - 96768 cos^° x + 130752 cos* x - 90048 cos^ x 



33264 cos^ X - 6300 cos^ x + 525) 



w. 



(3) 



154143 
22400 



(3.25) 
(3.26) 



Further comparing the above with (13. 6 p and (13. 8p . we identify the condensates 
perturbatively in Pq = A (see Table [1]). 



(a;(^V4'¥ 




Figure 1: (Colour online) Condensate values fl3.6p and f l3.8p for level j = 
{0,1,2,3} ( {purple,green,blue, orange} curves) boson stars. The red lines 
represent perturbative in p^ approximations, see Table [TJ 



3.3 Numerical boson stars 

In the previous section we identified the first four levels of a boson star 
perturbatively in the amplitude Pq. Here we report results for {u, (io,p|, d'^}, 
as well as {M, Q}, for generic p^ from the numerical solution of (13 .4^ . These 
results are collected in Figs. [H and El 

The top-left panel of Fig. [1] presents the frequency uj^^^ of a level-j boson 
star rescaled to that of a level-j oscillon frequency Uq (see (I3.13P ). We use 
purple/green/blue/orange color coding to denote j = 0---3. As pg (and 
correspondingly the mass and the charge — see Fig. [2] ) of a boson star 
grows, its frequency decreases. The remaining panels in Fig. [1] present the 
dependence of {dQ,pl,dl} as a function of pg. Notice that p| saturates; this 
saturation is the main obstacle in generating boson stars with ever increasing 
values of Pq (or mass). The red curves indicate perturbative approximations 
in Pq as collected in Tabled! (Fig. [1]), and perturbative approximation (13.141) 
in Q (right panel of Fig. [2]). 

In the limit of vanishing charge Q, the level-j boson star radial profile 
(j)^^'> is a single level-j oscillon (see (I3.13p . (I3.14p ) : 



(j)^^\x) (X ^/Q ej{x) . (3.27) 

For finite Q all the oscillons are excited. In Fig. |3] we present the spectral 
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Figure 2: (Colour online) Mass vs. Pq (left panel) and vs. charge (right 
panel) for level j = {0,1,2,3} ( {purple, green, blue, orange} curves) boson 
stars. The red lines represent perturbative in Q approximations, see (I3.14p . 

decomposition in the oscillon basis of the most massive level-j boson stars 
that we were able to construct 



.(i) 



7r/2 



dx (l)^^'{x)ei{x) tan^ x 



(3.28) 



SJ) 



Note that the maxima of c[ are achieved for i = j, much like in the small-Q 
limit. For all levels considered q"' approach a universal fall-off: 



cj^'^ oc (l + ^)-^ i>j 
represented by a dashed black curve in Fig. |3l 



(3.29) 



3.4 Perturbative stability of boson stars 

In this section we explore the linearized stability of boson stars. Consider 
perturbations of stationary solutions 03. 3p to leading order in e: 



0i(x,t) +i(t)2{x,t) = COS ^x( 0(x) + e(/i(t,x) - i(t){x)gi{t,x)) ]e 

A(t, x) = a{x) + e ai{t, x) , 
6{t, x) = d{x) + e 6i{t, x) . 



Further introducing 

/i(t,x) =Fi(x)cos(xt), 



Auit 



(3.30) 

(3.31) 



gi{t,x) = -Gi{x) sin(xt), 
the equations for ai(t,x) and 6i{t,x) can be solved explicitly: 

ai{t,x) = sin(2x)a(x) (w 0(x)^ G[{x) — (f)'{x) Fi(x)) cos{xt) , (3.32) 
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Figure 3: (Colour online) Spectral decomposition of level j = {0,1,2,3} 
({purple,green,blue, orange}) boson stars in oscillon basis, see f l3.29p . For 
large i all the curves approach a universal fall-off cf oc (1 + i)~^ (the black 
dashed curve). 



di{t,x) 



-2d{x) 



a{xY cos(x)0(x)G"/(x) sin(x) 



cos{x) (l){x)ijj sin(a:) 

- 2 cos(x)^0(x)^e^'^^^'^ sin(x)^a;^ + a{x){2a{x) cos(x)0'(x) sin(x) 
+ 2a(x) cos(x)^0(x) — a(x)0(x) — 2cos(x)^0(x) + 30(x))) G\{x) 
+ cos(x)0(x)e^'^^''^ G'i(x)x^sin(x) - 2cos(x)F(x)(cos(x)^0(x)0'(x) 

- cos(x)0(x)0'(x) - sin(x))a;e^'^^'''^ j cos(xt) • 

where Fi (x) and Gi (x) satisfy a coupled system of equations 



= F'; + -' ^"^^^^' + ' ^;^^^,^'" + ^ - V- - 20c. G'[ 



2uj 



(3.33) 



asm X cos X 



asin(x) cos(x) 



2 sinfx){/)^(/)' cos(x)'^ 



3sin('x)(/)^(/)'cos('x) 



2cos(x) 



a0 — 



2a(/)cos(x)^ — 2sinfx)a(/)' cosfx) 



G\ 



4acos(x)2(0')'-e'V 



6(0')'a + 3e^'^cu^ a-^Fi 



(3.34) 
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= G'l + I —6 cos(x)^0 + Qa(j) cos(x)^ + 2 sin(x)a0' cos(x) + 90 



- 7a(f) I ( a(f) sin(a;) cos(x) j G" + f 5a-^(f) + Aa'^ sin{x)(f)'(j) cos(x)'^ 

- lOa^ sin(x)0'0cos(x) — Aa^ sm{x)(f)' (p cos{xY + Ga^ sm{x)(f)' (p cos{x) 
+ 2e'^'^a(f)'^uj'^ cos(x)^ - e^'^acos(x)V^X^ - 2e'^'^a(j)'^uj^ cos(x)2 

+ e^'^a cos(x) V'x' + 2a\(f)'y cos(x) V' - 5a2(0')2 cos(x)V' 
+ 3a\(f)'y cos(x)202 + ge^^^ cos{xy(j)^uj^ + Qe^'' cos(x) ^V 

- ISe^'^ cos(x) V^w^ - 16a^02 cos(x)2 + 8a^0^ cos(x)^ + 9a0^ - 12aV^ 
+ 24a2 cos(x) V^ - 12a^ cos(x) V^ - 12a cos(x) ^^ + Aa cos(x)^02 

+ 2a^(0')^ cos(x)^ - 2a^(0')^ cos(x)M ( a^^^ sin(x)^ cos(x)M G'l 



+ ^£!!^ F{ + (402 cos(x)2 - 2a - 6</.2)e2wf a'^'l i^i 
a^0 V / 



(3.35) 



Notice that (I3.34p - p.35p are left invariant under the shift 

Gi^Gi + g, (3.36) 

where Q is an arbitrary constant. From (13.331) is it clear that this constant 
is fixed uniquely requiring that 

lim 5i(t,x)=0, (3.37) 

x— >-7r/2 

i.e., we keep the time coordinate at the boundary fixed. We do not have 
to worry about the shift symmetry (I3.36p . provided we rewrite (I3.34p - (l3.35p 
using 

dGi{x) = G\{x) . (3.38) 

Eqs. ( I3.34p -( j3735l) must be solved subject to constraints that Fi(x) and 
Gi{x)(l){x) are regular for x e [0,7r/2), and have only normalizable modes as 

X — 7- 7r/2, i.e., 

Fi (X p^ , Gi (X const , as p — )■ . (3.39) 

The latter regularity condition implies that Gi can have a simple pole (or 
dGi can have a double pole) precisely where (j){x) has a zerco. These poles 
represent a technical difficulty in identifying the fluctuations about excited 
boson stars — specifically, a straightforward shooting method: integrating 



''^Recall that excited level stationary boson stars are characterized by the number of 
nodes in the radial profile 4'{x), see p.3p . 
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from both boundaries with suitable boundary conditions and demanding con- 
tinuity at an arbitrary radial location will invariably encounter these poles 
rendering this method delicate to apply. 

Here we discuss the fluctuations about the lowest-level (ground state) bo- 
son stars, and also present analytic results for fluctuations about perturbatively- 
light excited boson starqj. 

Finally, since the system of equations fl3.34p -f l3l35|) is linear, we can fur- 
ther fix the normalizable mode Fi. 



lim ^ = 1 . (3.40) 



Fi 

Given (I3.39P and (I3.40p . we can specify the boundary conditions for 
{dGuF,}: 

■ at the origin of AdS, i.e., as x — )■ 0, 

dG^ =2g\x + 0{r') , F^ = f^ + 0{x^) ; (3.41) 

■ at the AdS boundary, z.e., as p — ?■ 

dGi =2gly + 0{x') , Fi = p^ + 0{p^) . (3.42) 

Note that along with x, the physical solution {dGi,Fi} is characterized 
by {Qii fQi92\ — which is the correct number of coefficients necessary to 
determine a unique (or isolated) solution for a pair of coupled second-order 
ODEs dSJl-driS]). 

Using the boundary conditions 03.411) and 03.42p it is easy to see that the 
charge of a fluctuating boson star does not change to leading order in e: 

.Qoc r^AP-"WW«WG;WUo. (3,43) 



jQ dx\ e°'(^) cos(a;)^ j 

3.4.1 Linearized fluctuations about light boson stars 

We report here the results for solving (13. 34p - ( 137351) for light boson stars, i.e., 
perturbatively in A, see section 13. 2[ In general, we search solution to above 
equations as a series 

dG['\x) = X dG'fl + A^ dG['l + C(A^) , 

Fi'^ = fS + A^ Fg + A^ Fi^l + 0{\^) , (3.44) 

X^'^ = xH^ + A^ x?^ + A^ xi'^ + 0(A^) ; 



^We verified explicitly that while dG\ to order 0{X^) has a double pole at the location 
of the zero of t/)^^^ (constructed perturbatively in A to order O(A^) inclusive — see P-lOp ). 
the full radial profile of physical fluctuations, G{ 4'^^\ is smooth for x € [0,7r/2]. 
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where j is the excitation level of a boson star. 

We finci^l: 
■ j = level, 



135 (0) 1215 2 


113892831 


(3.45) 


32 ' ^^ 128 ^ 


1254400 ' 


F{ Q = cos^ X , 




(3.46) 



.(0) cosix 



i3 



Fl"> = — — ^^--(5760cos(a;)^sin(x) - 14896 cos(x)^ sin(x) 
4480 sin(a;) 

+ 18738 cos(x)^ sin(x) + 7560 cos(a;)^x — 9153 cos(a;)^ sin(x) 

- 15120a: cos(x) + 1890sin(x)7r2 - 7560^^ sin (x)) - 2cos(x)^(32cos(x)^ 

-54cos(x)2 + 27)Ci,i, 

(3.47) 

dcfl = - —— sin(x) cos(2:)(5760 cos(a:)^ - 9376 cos(x)^ + 9900 cos(x)2 
' 4480 

- 2799) + 2 sin(x) cos(x) (8 cos(x)^ - 3) (4 cos(x)^ - 3)Ci,i . 

(3.48) 

The integration constant Ci.i is not fixed at order (9(A^), but is uniquely 
determineqil at order 0{X'^): 

Cu = i^. (3.49) 



j = 1 level, 



j = 2 level. 



j = 3 level. 



421120 



Xi^^ = 10, X« = -^. (3.50) 



X?^ = 14, ^f) = J_^^, (3.51) 



(3) 10 (3) 62451 

X|,^ = 18, X\' = -^^. (3.52) 



3.5 Linearized fluctuations about j = boson stars 

The spectrum of linearized fiuctuations about j = boson stars is presented 
in Fig. m We find that over the whole range of charges Q we were able 



^For excited levels we present only the coefficients {x^"' } , i = 0, 1 , j = 1 • • • 3. 
'^This pattern extends to higher orders in A: a solution at order A^", {i^i,2n , Gi^2n-i}, 
is determined up to a constant Ci_„, which is being uniquely fixed at order ©(A^'^""'"^)). 
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Figure 4: (Colour online) Spectrum of linearized fluctuations about j = 
boson stars as a function of pg (black dots). The solid orange/red curves are 
successive approximations to x^ = (x(Po — '^)) i^ -^^^ see (I3.44p . 

to construct j = boson stars, the frequency of their fluctuations squared 
(x^) is positive. This strongly suggests that the ground state boson stars are 
perturbatively stable. 

As discussed in the previous section (see (I3.50|) - O3.52|) ). excited boson 
stars are perturbatively stable for small charge. Our numerical simulations 
suggest that both the j = and the excited boson stars are nonlinearly 
stable. 

4 Nonlinear Results 

We take the constructed boson star solutions described above {(f){x), d{x), a{x),u} 
and employ them to provide initial data for our dynamical studies via 

5-, (4.1) 



cos^x 



^^ = ^^Sl, (4.2) 

cosx 



u0e^ 2 
a 



n, = ^^Sl (4.3) 



where the metric functions are obtained by solving the constraints. We con- 
firm convergence of the obtained solutions (by monitoring the constraints 
residuals, charge and mass conservation and self-convergence vs time) as res- 
olution is increased (see also [T3]). 

16 



Perturbed, Genuine Boson Stars: We concentrate on studying the 
behavior of these solutions when perturbed, and have considered various 
forms of perturbation with quahtatively similar results. For concreteness, we 
here present results obtained with Gaussian perturbations parametrized as 

G{x) = ee"*-'""^"^ '^ and add it to the boson star solution via 

0. = [(j)/cos^x + G{x)]6l, (4.4) 

$. = [(f)' / cos x + G'{x)]6l, (4.5) 

n, = [u(l){e''/a)+G']5l (4.6) 

In analogy with previous studies, we set the amplitude of the Gaussian per- 
turbation to e. We note that because the constraints are solved numerically 
at the initial time for a{x, 0) and S{x, 0), this perturbed initial data together 
with the obtained metric variables is fully consistent with the equations of 
motion. Because our numerical implementation is fully non-linear, it natu- 
rally probes the nonlinear behavior of this system, and we are particularly 
interested in the stability of these boson star solutions. 

Recall that Ref. [1] found that pulses of scalar field in AdS are unstable 
to black hole formation, and subsequent arguments in [5] supported this view 
of generic instability. A reasonable expectation in light of those works is that 
any perturbation of the boson star will behave in a similar way; that is, one 
expects such a perturbation to travel back and forth between origin and AdS 
boundary, sharpening with each pass, leading eventually to BH formation. 
However, our early studies suggested quite the opposite (for sufficiently small 
perturbations); that instead boson stars were stable [17] prompting us to 
study this system more broadly and deeply. 

For small perturbations, very long-lived, regular solutions were obtained 
describing perturbed boson stars. For such long-lived solutions, we monitored 
the metric functions (e.g. max|l — A{x,t)\ and 6{0,t)) and they showed 
no signs of instability to BH formation. These observations prompted a 
thorough study of the instability in AdS; independent work via perturbative 
studies also pointed out that AdS should be stable for several families of 
solutions [8]. 

A variety of boson star solutions were studied, including members of 
levels 0, 1, and 3 (level 2 solutions presented regularity issues near the AdS 
boundary and we defer such analysis for future work). All examples appeared 
stable. Interestingly, in asymptotically flat scenarios, excited boson stars 
are generally unstable, radiating energy and settling into a ground state 
solution [20]. In AdS, however, there is no way to rid itself of excess charge, 
which presumably explains their stability. Nevertheless, this property of AdS 
does not explain how the boson star can be immune to the weakly turbulent 
instability. Below, we present an argument to this end, but flrst we discuss 
the behavior of a different family that lends support to our argument. 
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Fake Boson Stars: Of course numerical evolutions are limited to finite 
times, and so one cannot rule out that instability will manifest after the code 
has been stopped or beyond the time for which one trusts the results. To 
better assess the observed behavior, we compare these long-lived solutions 
to a different family which can be considered "nearby" in some sense. This 
family, which we refer to as fake boson stars, represents purely real initial 
data with the same mass and profile as their counterpart genuine boson stars. 
A fake counterpart of some boson star solution of (14. 3 p is achieved via the 
transformation 

^fake^^BS^ j^fake^j^BS^ 0^^'^'= = . (4.7) 

Remarkably, the evolution of this family also yields regular, long-lived 
solutions for small perturbations that do not collapse to a black hole. Figs. [5] 
and [6] illustrate the time of collapse as a function of e both for genuine and 
fake boson stars. As indicated in the figures, successively higher resolutions 
largely coincide with differences only apparent at the latest times. In all 
cases, the results indicate collapse times increasing quickly as e decreases 
with no signs of collapse for smaller amplitudes of perturbation. 

Notice that these fake solutions are not stationary and have no charge, 
two seemingly essential features of genuine boson stars, and so their apparent 
immunity to this weakly turbulent instability is surprising. This "stability" 
is apparently not tied to special features (e.g. charge or stationarity) but 
instead suggests that the dynamics undergoes something akin to a frustrated 
resonance in which amplitudes increase at times but then disperse. In par- 
ticular, one essential aspect common to both genuine and fake boson stars 
appears to be their non-compact, long-wavelength nature. Because they have 
energy distributed throughout the domain, modes no longer propagate co- 
herently. Instead there is a continuing competition between dispersion and 
gravitational contraction; collapse to a black hole or not is then determined 
by the outcome of this competition. 

Large a: Admittedly, this argument is far from rigorous. But if it holds, 
then it would imply many other forms of stable initial data. In particular, 
perhaps other forms of initial data may be immune to this weakly-turbulent 
instability when its extent is large. To explore this conjecture, we adopt the 
same form of data considered in many previous studies of this instability (such 
as those in [H [131 [IHl Ell [22]). We thus consider this family again, which 
takes the following form in our rescaled variables 

$i(0,a:) = 0, Ui{0,x) = —e-^^^^ cos^-'^x Si. (4.8) 

vr 

To test the possibility of regular development of this data, we considered 
the time development of Eq. (14.81) with varying values of a. The results are 
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Figure 5: (Colour online) Collapse times for Gaussian perturbations of a 
ground state boson star (0i(O,O) = 0.253) and its corresponding fake star. 
Increasing resolutions are shown. For short collapse times, resolutions agree. 
However, for the longest evolutions, higher resolutions are needed. Even with 
very high resolutions, small e evolutions show no sign of collapse. 
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Figure 6: (Colour online) Collapse times for Gaussian perturbations of a first 
excited state boson star (0i(O,O) = —0.272) and its corresponding fake star. 
As in Fig. |5l higher resolutions are also shown with differences among the 
resolution appearing only at very late times. 
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plotted in Figs. [7] and [8] which show the time of collapse as a function of e 
for various values of a. As is evident from the figures, for small values of a 
[a < 0.3) the collapse time increases monotonically as e decreases; however 
for larger values of a the collapse time increases abruptly as e is decreased. 
Notice that this abrupt growth in collapse time behaves quite similarly to 
that seen for boson stars and fake stars, suggesting that for sufficiently large 
(T, the behavior would be regular. Furthermore, an analysis of the Fourier 
power spectra of cases below a ~ 0.3 reveal that the spectra monotonically 
shift to higher frequencies as time progresses. In contrast, for cases above 
a K. 0.4 they do not do so. Instead the shift saturates and the spectral 
content oscillates within a narrow window of frequencies. 

That large-(T initial data is immune to the weakly turbulent instability 
is consistent with the argument that widely distributed mass energy pre- 
vents the coherent amplification typical of the instability. It is interesting 
to consider what would happen in the semilinear wave equation on a fixed 
AdS background as studied in [22]. That model shows many of the same 
characteristics as the gravitating scalar collapse, but the nonlinear potential 
plays the role of the attractive, focusing effect that gravity plays here. How- 
ever, numerical evidence from that model suggests that there is no large-a 
effect, lending support to the idea that the distributed mass-energy affects 
the spacetime in a way that disturbs the coherent amplification. 

A change in behavior such as this merits a closer examination of the 
"transition region" between apparent stability and black hole collapse. Fig. [9] 
illustrates this region 0.3 < o" < 0.4 in more detail. Interestingly, in this tran- 
sition region the time-to-collapse exhibits a seemingly oscillatory behavior 
prior to displaying the characteristic rapid growth as e is decreased. 

From Figs. [7]and[9l it is clear that for a > 0.4 there is some emin below 
which initial data does not form a black hole. The idea of this function 
emin(o") is similar to emin(2;max) studied in [22]. Preliminary study of eram{cr) 
shows it to be a roughly exponentially decreasing function (after the apparent 
discontinuity at a ~ 0.4 in which emin = — )■ 5.6). The behavior below emm 
is demonstrated in Fig. [Ill In particular, for this weak initial data when the 
metric is frozen at its initial profile, the evolution demonstrates dispersion. 

It is instructive to study the spectral decomposition of the initial data 
fl4.8j) in the oscillon basis for different a. To relate with the analysis in Fig. [TJ 
we keep ae = 1 fixed. For a select set, i.e., a = {0.0625,0.1,0.2, ■ ■ -0.7}, we 
compute the spectral coefficients Ci(cr), see (12.51) . 



dx tan^x A{0,x)e~^''^'''^Ui{0,x) ei{x) 





u(^ 



(4.9) 



where A{0,x) and 6{0,x) are obtained from integrating (12. 8p with initial 
data (14. Sp . The resulting spectral decompositions are collected in Fig. ^U\ 
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Comparing with the spectral decomposition of boson stars (see Fig. |3]), here, 
the large-j decay of the spectra is approximately exponential, instead of a 
power-law as in (I3.29p . The spectral profile achieves a minimum around 
a ~ 0.2 — 0.4, which is roughly the critical value of a separating the stable 
and unstable regions in the parameter space of the initial data (14. 8p . 
Restricted domain: The last nonlinear effect presented in this section 
concerns evolutions conducted in a restricted domain. In particular, an artifi- 
cial, reflecting boundary condition is applied at some x^ax < 7i"/2, restricting 
the propagating pulse to some subdomain of AdS. The motivation for this is 
to study whether this turbulent instability is itself just the manifestation of 
the nonlinear attraction of gravity occurring in a bounded domain, or instead 
some particular property of the full AdS (and hence would be destroyed by 
this restriction). 

As found in [13], the imposition of such a reflecting boundary condition 
does not eliminate black hole formation after multiple bounces. However, a 
minimum value of e was found, below which no such black hole formation 
occurred. In the semilinear model, it was found that the boundary condition 
resulted in dispersion not seen with the full AdS domain. 

Here, we revisit this problem, showing the customary time of collapse plot 
in Fig. [121 As mentioned, there is some similarity with the large-a effect in 
the existence of some emin- However, in contrast, one does observe the "stair- 
step" decrease in collapse time for e increasing above emm, characteristic of 
successive bounces. 

It is also possible that such restricted-domain evolutions result from two 
effects: (i) the imposition of the reflecting wall introduces dispersion as in the 
semilinear model, and (ii) as the domain shrinks, there may be some effect 
due to the fact that, for fixed a, the fractional support of the initial data is 
increasing. 

5 Conclusions 

We have constructed boson star solutions in global AdS and shown that 
they are stable at linear order. Numerical studies of their dynamics strongly 
suggest that these boson stars, both ground state and the first few excited 
states, are non-linearly stable. Along with solutions presented in [15], there is 
now considerable evidence that the instability of AdS to scalar perturbations 
reported in [1] is limited in scope; that there exist non-trivial, dynamical 
examples of stable solutions in AdS. These results are consistent with the 
perturbative arguments of |H] for the stability of boson stars, geons and 
solitons. 

Comparison of the lifetimes of perturbed boson stars with other, non- 
stationary solutions, our fake boson stars, reveal an even wider class of ini- 
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Figure 7: (Colour online) Collapse times for initial data of the form Eq. (14 .Sp 
with varying width values, a. Because changes to a affect the amount of 
mass, the natural parameter against which to plot is ae not just e (also see 
Fig. IH] for this data plotted versus e). Note that for a < 0.3 the standard 
behavior is observed where collapse eventually occurs for any e. In contrast 
for a > 0.3, there appears to exist a threshold e* below which collapse does 
not occur. For initial data above the transition, a > 0.3, evolutions with 
smaller e than shown reached at least t ^ 2000 with no signs of eventually 
collapse. 
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Figure 8: (Colour online) Same data as Fig. [7] but with an un-rescaled ab- 
scissa. As one's eye moves from upper-right down to lower-left, the initial 
data parameter a increases and the behavior of the time of collapse, tc, 
changes dramatically from the "usual" stair-step to something else entirely 
with a very sharp transition. Note that the spacing between runs is very 
non-uniform and that tuning is necessary to resolve the transition. 
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Figure 9: (Colour online) Collapse times for a range of a that demonstrates 
the transition from turbulent BH formation to frustrated resonance. For 
a = 0.3, collapse appears inevitable for any value of e in contrast to the 
results for a = 0.4. Interestingly, a "bump" appears for these values of a 
in which the collapse times demonstrate a lack of monotonicity. Beginning 
with a = 0.3 (magenta, solid squares) around e ~ 8, one sees a small bump 
that, as one looks to higher a, sharpens and occurs at smaller e values. 
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Figure 10: (Colour online) Spectral decomposition of the ini- 
tial data gSD with a = {0.0625,0.1,0.2,0.3,0.4,0.5,0.6,0.7} (yel- 
low,blue,cyan,niagenta,black,red, green, orange) and ere = 1 in the oscillon 
basis. Spectral coefficients Ci{a) (see (14. 9p ) start relatively large at high os- 
cillon numbers for a = ^ (yellow); they decrease for o" = 0.1 (blue), and 
achieve a minimum profile around a = 0.2 (cyan) or cr = 0.3 (magenta); then 
they increase for a = 0.4 (black), a = 0.5 (red), a = 0.6 (green) and a = 0.7 
(orange). The minimum of the spectral coefficients proffie roughly coincides 
with the critical value of cr ~ 0.2 — 0.4 which separates stable and unstable 
regions in the parameter space of the initial data (14. 8p . 



tial data which appears to frustrate the resonance and thereby avoid collapse 
for sufficiently small amplitude. Indeed, using the very same family of the 
seminal work of [1], numerical evolutions suggest that for initial data with 
a ^ 0.4 the instability can be avoided. More generally, initial data with 
widely distributed mass-energy appears to be similarly immune to the turbu- 
lent instability. This behavior is in contrast to what is observed in solutions 
to the semilinear wave equation in (flat) AdS. There, a singularity always 
form [22], even for large-o", indicating that gravity plays a key role in the 
dynamics. In particular, a heuristic argument suggests that the widely dis- 
tributed mass-energy distorts the space sufficiently to introduce dispersion 
and thereby oppose the concentrating effect of the instability. 

The picture that emerges is a phase space with (at least) two regions. 
One region is subject to the weakly turbulent instability and therefore col- 
lapses for any initial "amplitude." The other region can avoid the turbulent 
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Figure 11: (Colour online) Demonstration of dispersion introduced by widely 
distributed mass-energy Shown is the behavior of 11 at the origin during the 
evolution of e = 1, a = 1 initial data. By keeping A = 1 and 5 = 0, we evolve 
the scalar field in a pure AdS background, and this results in a periodic 
solution (dashed, red line). Instead, by solving for the initial metric A{x,0) 
and 6{x,0), we evolve the scalar field on a fixed background (solid, blue 
line). This fixed-background evolution displays the dispersion introduced by 
the fixed metric. The fully dynamical evolution (dashed magenta) is also 
shown. In the intervening time (see inset), this solution shows the periods of 
focussing and dispersal typical of what we call frustrated resonance. Finally, 
at late times, we copy the pure AdS solution from the period < t < vr 
and display it shifted in time (dashed solid); that it overlays the pure AdS 
solution shows that scalar solutions on a background of AdS are periodic. 
This figure is similar to Fig. 5 of [22] which shows the dispersion introduced 

by a restricted domain. 
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Figure 12: (Colour online) Collapse times of the initial data in Eq. f l4.8p with 
a = 1/16 with an artificial, reflecting wall at a various positions, Xmax- Also 
shown (solid black) are the results for the full domain with no reflecting wall. 
Note that the finite domain evolutions demonstrate a threshold emin below 
which collapse does not occur. 



28 



instability and contains oscillons, boson stars, geons, and similar symmetric 
solutions. However, it appears this second "stable" region contains a wider 
class of solutions with no periodicity or stationarity, namely fake boson stars, 
large-o" and similarly distributed families of initial data. To be clear, this sec- 
ond region need not be strictly stable and can certainly possess collapsing 
solutions. What is important is that in this second region, one can choose 
a sufficiently small amplitude such that the weakly turbulent instability is 
avoided. 

We do not know at this stage what is the precise criterion that separates 
the parameter region of scalar field initial data resulting in BH formation 
from the region of non-linear stability. It can be argued that such a criterion 
is encoded in the spectral decomposition of the initial data in the oscillon 
basis. Indeed, consider the region (assuming it existqfl) where collapse occurs 
for arbitrarily small amplitude of the scalar field. In this limit, the full initial 
data is the oscillon spectrum, as the backreaction can be safely ignored. In 
this paper we presented strong evidence for initial configurations that do not 
collapse in the limit of vanishingly small amplitude. Thus, the distinction 
between stable and non-stable configurations (at least for small amplitudes) 
must be hidden in the initial scalar field spectral data. We have seen that 
initial profiles of [1] (14. 8 p become stable for small e as a increases; the latter 
increase results in softening the decay of the asymptotic oscillon spectral 
coefficients (see Fig. [TOj) . Likewise, boson stars have an asymptotic power-law 
oscillon spectral decomposition, in contrast to the exponential-decay profile 
for initial data (gSD (see Fig. [3]). 

It is important to further investigate the nonlinear stability of AdS. The 
issue has profound implications for a dual boundary conformal field theory, 
as it identifies CFT initial configurations that fail to thermalize. There is 
no obvious symmetry criterion "protecting" such configurations. A possible 
future direction is to investigate the collapse of initial configurations specified 
by their oscillon spectral decompositions with various trial profiles Cj. The 
primary goal, of course, is the identification of the stability criteria — it 
is possible that the latter can be established from analysis of the weakly 
nonlinear regime only. Finally, it is interesting to analyze the stability of 
more general configurations — such as boson stars with local bulk charge, as 
recently discussed in 
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